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Abstract. 

We present a reduction procedure for gauge theories based on quotienting out the 
kernel of the presymplectic form in configuration-velocity space. Local expressions 
for a basis of this kernel are obtained using phase space procedures; the obstructions 
to the formulation of the dynamics in the reduced phase space are identified and 
circumvented. We show that this reduction procedure is equivalent to the standard 
Dirac method as long as the Dirac conjecture holds: that the Dirac Hamiltonian, 
containing the primary first class constraints, with their Lagrange multipliers, can be 
enlarged to an extended Dirac Hamiltonian which includes all first class constraints 
without any change of the dynamics. The quotienting procedure is always equivalent 
to the extended Dirac theory, even when it differs from the standard Dirac theory. 
The differences occur when there are ineffective constraints, and in these situations we 
conclude that the standard Dirac method is preferable — at least for classical theories. 
An example is given to illustrate these features, as well as the possibility of having 
phase space formulations with an odd number of physical degrees of freedom. 
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1. Introduction 

The dynamics of gauge theories is a very wide area of research because many 
fundamental physical theories are gauge theories. The basic ingredients are the 
variational principle, which derives the dynamics out of variations of an action 
functional, and the gauge principle, which is the driving principle for determining 
interactions based on a Lie group of internal symmetries. The gauge freedom exhibited 
by the complete theory becomes a redundancy in the physical description. The classical 
treatment of these systems was pioneered by Dirac (1950, 1964), Bergmann (1949), 
and Bergmann and Goldberg (1955). Dirac's classical treatment in phase space (the 
cotangent bundle for configuration space) has been shown (Gotay and Nester 1979, 
1980, Batlle et al 1986) to be equivalent to the Lagrangian formulation in configuration- 
velocity space (the tangent bundle). One ends up with a constrained dynamics with 
some gauge degrees of freedom. One may choose, as is customary in many approaches 
(Pons and Shepley 1995) to introduce new constraints in the formalism to eliminate 
these unwanted — spurious — degrees of freedom. This is the gauge fixing procedure. 

There are approaches other than gauge fixing. For instance, the method of Faddeev 
and Jackiw (1993)and Jackiw (1995) is to attempt to reduce the system to its physical 
degrees of freedom by a process of directly substituting the constraints into the canonical 
Lagrangian. It has been proved (Garcia and Pons 1997) that, as long as ineffective 
constraints — functions that vanish in the constraint surface and whose differentials also 
vanish there — are not present, the Faddeev- Jackiw method is equivalent to Dirac's. 

A reduction procedure (Abraham and Marsden 1978, Sniatycki 1974, Lee and Wald 
1990) which seems to be particularly appealing from a geometric point of view consists 
in quotienting out the kernel of the presymplectic form in configuration- velocity space in 
order to get a reduced space, the putative physical space, with a deterministic dynamics 
in it, that is, without gauge degrees of freedom. One must be careful that these 
techniques do not change the physics, for example by dropping degrees of freedom, 
and that they are applicable in all situations of physical interest. For example, we know 
of no treatment of this technique which applies to the important case when there are 
secondary constraints — one purpose of this paper is to provide this treatment. 

In this paper we present a complete reduction method based on quotienting out 
the kernel of the presymplectic form. We establish a systematic algorithm and prove 
its equivalence with Dirac's method, but only so long as ineffective constraints do not 
appear. Our procedure turns out to be equivalent to Dirac's extended method, which 
enlarges the Hamiltonian by including all first class constraints. It differs from the 
ordinary Dirac method (supplemented by gauge fixing) when ineffective constraints 
occur. Since the ordinary Dirac method is equivalent to the Lagrangian formalism, it is 
to be preferred in classical models. 

We will consider Lagrangians with gauge freedom. Thus they must be singular: 
The Hessian matrix of the Lagrangian, consisting of its second partial derivatives with 
respect to the velocities, is singular, or equivalently, the Legendre transformation from 
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configuration-velocity space to phase space is not locally invertible. Singular also 
means that the pullback under this map of the canonical form u> from phase space 
to configuration- velocity space is singular. 

In order to proceed, we first compute, in section 0, in a local coordinate system, a 
basis for the kernel of the presymplectic form. Our results will be in general local; global 
results could be obtained by assuming the Lagrangian to be almost regular (Gotay and 
Nester 1980). In section ^ we will single out the possible problems in formulating the 
dynamics in the reduced space obtained by quotienting out this kernel. In section [| we 
will solve these problems and will compare our results with the classical Dirac method. 
It proves helpful to work in phase space here, and we end up with a reduced phase space 
complete with a well-defined symplectic form. In section ||we illustrate our method with 
a simple example (which contains ineffective constraints). We draw our conclusions in 
section ||. 

2. The kernel of the presymplectic form 

We start with a singular Lagrangian L(q\ q % ) (i — 1, • • • , N). The functions 
pi(q,q) ■= dL/dq 1 

are used to define the Hessian W%j = dpi/dqi, a singular matrix that we assume has 
a constant rank N — P. The Legendre map TL from configuration-velocity space (the 
tangent bundle) TQ to phase space T*Q, defined by pi = pi(q,q), defines a constraint 
surface of dimension 2N — P. 

The initial formulation of the dynamics in TQ uses the Lagrangian energy 

E L := p4 - L 

and X, the dynamical vector field on TQ, 

zx^l = d(E L ) , (2.1) 

where 

u; L := dq s A dp s 

is the pullback under the Legendre map of the canonical form u> = dq s A dp s in phase 
space, cjl is a degenerate, closed two-form, which is called the presymplectic form on 
TQ. In fact there is an infinite number of solutions for equation ( [2.1D if the theory has 
gauge freedom, but they do not necessarily exist everywhere (if there are Lagrangian 
constraints). X must obey the second order condition for a function: 

Xg* = ?' ^X = g*A + ^M)A, 

where A s is partially determined by equation Q2.1| ). 

At first sight, the kernel of wl describes, in principle, the arbitrariness in the 
solutions X of equation fl2.1p . Therefore it is tempting to think that in order to 
construct a physical phase space, we must just quotient out this kernel. The complete 
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implementation of this procedure, which we are presenting in this paper is, first, far 
from obvious, and second, as we will see, fraught with danger. 
Let us first determine a basis for 

/C := Ker(u;L) 
in local coordinates. We look for vectors Y satisfying 

i Y u h = . (2.2) 
With the notation 



Qqk Qgk 

equation (|2.2j) implies 

eWy = , (2.3a) 

e i A ij -p i W ij = , (2.36) 

where 

A - d& d Pj 
13 ' dqi dq* • 

Since W is singular (this causes the degeneracy of it possesses null vectors. It 
is very advantageous to this end to use information from phase space. It is convenient 
to use a basis for these null vectors, 7*, (/x = 1, . . . ,P), which is provided from the 
knowledge of the P primary Hamiltonian constraints of the theory, (f>y. Actually, one 
can take (Batlle et al 1986), 

where TL* stands for the pullback of the Legendre map JFL : TQ — > T*Q. According 
to equation ( |2.3q| ), e l will be a combination of these null vectors: e* = A M 7^. Notice that 
we presume that these primary constraints are chosen to be effective. 

To have a solution for (3 l we need, after contraction of equation ( |2.36| ) with the null 
vectors 7^, 




= A'T^Vrf = aw HM hrf - lS ^ 



dqi dq % J \ dpj 
which is to be understood as an equation for the A^s. We then use the identity 

, (2.5) 




which stems from the fact that (q, p) vanishes identically; we get 



(i)^(i) ^W^W 



= 



dpi dq 1 dqi dpj 
AW «/>«}) . (2.6) 
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Condition (|2.6| ) means that the combination X^ip^ must be first class. Let us split the 
primary constraints </>]L between first class <f/^ and second class <f)^) at the primary 
level, and we presume that second class constraints are second class everywhere on the 
constraint surface (more constraints may become second class if we include secondary, 
tertiary, etc, constraints). They satisfy 

{0£?,0i 1) }=pc, det|{0« 0«}|^O, (2.7) 

where pc stands for a generic linear combination of primary constraints. Equations ( |2.6| ) 
simply enforce 

= . 

Consequently a basis for the e l will be spanned by 7 W , so that 
= ^% 

for A Ml arbitrary. Once e l is given, solutions for ft will then be of the form 

ft = ^pu + , 

where the ^ are arbitrary functions on TQ. We will now determine the ft . 

To compute ft it is again very convenient to use Hamiltonian tools. Consider any 
canonical Hamiltonian H c (which is not unique), that is, one satisfying Ei, = J-"L*(H C ). 
Since we know from the classical Dirac analysis that the first class primary constraints 
4$} may produce secondary constraints, 

*g? = {*£Uc}, 

we compute (having in mind equation fl2.3£| )) 

< A , g^g) 



13 

-• S I g(^) (Of* 

%Ai + > (2-8) 



where we have used the operator if defined (Bathe et al 1986, Gracia and Pons 1989) 
by 

\ oq l J oq l \ dpi 
This operator satisfies (Batlle et al 1986, Pons 1988) 

Kf = TV ({/, H c }) + v*{q, q)TL* ({/, $>}) , (2.9) 
where the functions are defined through the identities 

<f = JL* ({g*, # c }) + ^(g, g) JL* ({g\ . (2.10) 

Property ( |2.9D implies, for our first class constraints, 
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which has been used in equation fl2.8|) . Let us continue with equation 



dqj 




= WijK-^- , (2.11) 
opt 

where we have omitted some obvious steps to produce the final result. We can read off 
from this computation the solutions for equation 



pi = K^L ^jr L *(^S_\ . (2 . 12) 



/,] dpj \ dpj 



Therefore, a basis for K, is provided by: 

.A 



and 



^ : =<7w+#i^- (2-135) 



9 a, 9 
h J 

dqj 1 ^ dqj 

Vectors in equation ( |2.13q| ) form a basis for Ker(TjFL), where TTL is the tangent 
map of TL (also often denoted by TL^. They also span the vertical subspace of K,: 
Ker(TjFL) = Ver(/C). This is a well known result (Carinena et al 1988), but as far as 
we know equations ( [2.134 |2 . 13 6| ) are the first explicit local expression for /C itself. 

All other results (Carinena 1990), obtained on geometrical grounds, for /C are 
obvious once the basis for this kernel is displayed, as it is in equations ( |2.134 |2.13fe| ). 
For instance, it is clear that dimJC < 2 dim Ver(/C). Also, defining the vertical 
endomorphism 

we have S(/C) C Ver(/C). The case when 
S(/C) = Ver(/C) , 

corresponds, in the Hamiltonian picture, to the case when all primary constraints are 
first class (indices p = indices p\). These are the so-called Type II Lagrangians (Cantrjin 
et al 1986). S(/C) = corresponds to the case when all primary constraints are second 
class (indices p = indices p^). 

Equation (|2.13c| ) implies, for any function / on T*Q, 

r M (^L*(/)) = 0- (2-14) 
The corresponding equation for A Ml is: 

A Ml (JL*(/))=^L*({/,0«}) . (2.15) 
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Since we will make use of this property below, we now prove this result. The action of 



We use equations Q2.36Q , 




A w (^L*(/)) 



= ({/,*#}) • 

The commutation relations (Lie Brackets) for the vectors in equations ( |2.13q| , |2.13£| ) 
are readily obtained, and we present these new results here for the sake of completeness. 
We introduce the notation 

(commutation of first class constraints is also first class). We arrive at 

[r M ,r„] =0, (2.16a) 
[T MI A M1 ] =JV(Al w )r v , (2.166) 

[A m , A„J = JL* A Pl + ^ JL* (B^^, T p , (2.16c) 

where the ?A are defined in equation ( p.lOj ). Observe that the number of vectors in 
equations (|2.134 |2.13fe|) is even because = — is the number of second class 



primary constraints (at the primary level), which is even. 

Because the algebra of fC is closed, the action of K on TQ is an equivalence relation. 
We can form the quotient space TQ/JC and the projection 

7T : TQ — ► TQ/JC . 

TQ/JC is endowed with a symplectic form obtained by quotienting out the null vectors 
of Wl (that is, wl is projectable to TQ/fC). The space TQ/fC is not necessarily the 
final physical space, however, because we have not yet formulated the dynamics of the 
system: We now turn to the question of the projectability of the Lagrangian energy. 

3. Obstructions to the projectability of the Lagrangian energy 

In order to project the dynamical equation (|2.1|) to TQ/fC, we need E L to be projectable 
under tt. However, in order for Ei, to be projectable we must check whether it is constant 
on the orbits generated by IC as defined by the vector fields of equations ( [2.134 P-136|) . 
Indeed T^Ei) = 0, but from equation (|2.15| ), 

A m (£0 = -^*W?) > 
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where 

IlTL*(<j))Z>) ^Ofor some /xi, then 0^/ is a secondary Hamiltonian constraint. As a side 
remark, note that in this case jFL*(0/ tl ) is a primary Lagrangian constraint. In fact it 
can be written (Batlle et al 1986) as 

tv ($>) = [L]al , 

where [L]j is the Euler-Lagrange functional derivative of L. 

We see that there is an obstruction to the projectability of E L to TQ/JC as long as 
there exist secondary Hamiltonian constraints or equivalently if there exist Lagrangian 
constraints. 

One way to remove this problem (Ibort and Marin-Solano 1992), Ibort et al 1993) is 
to use the coisotropic embedding theorems (Gotay 1982), Gotay and Sniatycki 1981) and 
look for an extension of the tangent space possessing a regular Lagrangian that extends 
the singular one and leads to a consistent theory once the extra degrees of freedom 
are removed. This method is equivalent to Dirac's, but only if there are no secondary 
Hamiltonian constraints. However, if there are, which is precisely our case, the dynamics 
becomes modified and thus changes the original variational principle. Instead of using 
this technique we will try to preserve the dynamics. 



4. Physical space 

In the cases where secondary Hamiltonian constraints do exist (for instance, Yang-Mills 
and Einstein-Hilbert theories), we must find an alternative reduction of TQ in order to 
obtain the projectability of E^. 

The initial idea was to quotient out the orbits defined by equations ( [2.134 |2 . 1 3 fcj) . 



Since T^E^) = we can at least quotient out the orbits defined by equation (|2.13o|) . 
But this quotient space, TQ /Ker {TTL\ is already familiar to us: It is isomorphic to 
the surface M\ defined by the primary constraints in T*Q. In fact, if we define tti as the 
projection 

tti : TQ — ► TQ/Ker(TFL) , 
we have the decomposition of the Legendre map J-L — i\ o 7Tl, where 

r £Q_ _ T^y-) 

1 ' Ker(TjFL) V ' 



with 

,( 



\Ker(TjFL) J 



M 1 



Now we can take advantage of working in Mi C T*Q. Let us project our original 
structures on TQ into M\. Consider the vector fields A Ml . Equation ( 2.15Q tells us that 



the vector fields A Ml are projectable to M\ and that their projections are just { — , ($}} 
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In fact these vector fields {— , $1} are vector fields of T*Q, but they are tangent to Mi 
because <f$ are first class (among the primary constraints defining Mi). Incidentally, 
note that the vector fields {— , 4>^)} associated with the second class primary constraints 
in T*Q are not tangent to M\. 

Formulation in Mi of the dynamics corresponding to equation ( |2.1D uses the pre- 
symplectic form u)\ defined by u>i = z*u>, where u> is the canonical form in phase space, 
and the Hamiltonian Hi defined by Hi = i*H c , with H c such that J-L*(H C ) = E^. The 
dynamic equation in Mi will be: 

i Xl Wi = dHi . (4.1) 

The null vectors for u)i are { — ,4>$} (more specifically, their restriction to Mi). 
(This result is general: The kernel of the pullback of the symplectic form to a constraint 
surface in T*Q is locally spanned by the vectors associated, through the Poisson Bracket, 
with the first class constraints among the constraints which define the surface.) To 
project the dynamics of equation (fO]) to the quotient of Mi by the orbits defined by 

{-<$}: 

V, := — * ■ («> 
({-*$}) 

we need the projectability of Hi to this quotient manifold. To check this requirement 

it is better to work in T*Q. Then projectability of Hi to V\ is equivalent to requiring 

that {<f>$, H c }\ Mi =0. 

Here lies the obstruction we met in the previous section, for it is possible that 
, H c }\ Mi 7^ for some constraints <p^. Let us assume that this is the case. As we 

did before, we define 

These constraints may not be independent, some of them may vanish on Mi, and some 
previously first class constraints may become second class. Those that do not vanish 
are secondary constraints and allow us to define the new surface M 2 C Mi (we define 
the map i 2 : M 2 — ► Mi) by <$} = 0. 

We can now form the projection of H 2 := %%Hi to M 2 /({ — , }), but the projection 
of uj 2 '■= i\<jJi can be still degenerate in this quotient space, since uj 2 may have acquired 
new null vectors (and may have lost some of the old ones). In fact, once all constraints are 
expressed in effective form, Ker(o;2) is generated under the Poisson Bracket associated 
with uj by the subset of effective constraints that are first class with respect to the whole 
set of constraints defining M 2 . If there is a piece in this kernel that was not present in 
Ker(u>i), then new conditions for the projectability of H 2 will appear. 

The dynamic equation in M 2 is 

«x 2 ^2 = dH 2 . (4.3) 

(2) 

It is still convenient to work with structures defined in T*Q. Suppose that (jy^ is any 
secondary, first class, effective constraint in M 2 ; therefore {— , 0J* 3 } £ Ker(u2) but 
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(2) (2) 

{— j 0M2 } ^ Ker(cJi). The new project ability condition for H 2 induced by (p^ is 

{C#c}L 2 =o. 

This means that we might find new constraints if this condition is not satisfied. A new 
surface M 3 will appear, and a new kernel for a new u>3 should be quotiented out, and so 
on. We will not go further because we are just reproducing Dirac's algorithm in phase 
space (Dirac 1950, 1964, Bathe et al 1986, Gotay et al 1978). We do have a shift of 
language, however: What in Dirac's standard algorithm is regarded as a condition for 
the Hamiltonian vector field to be tangent to the constraint surface is here regarded as 
a projectability condition for the Hamiltonian to a quotient space. 

To summarize: The constraint surface Mi is defined by the primary constraints 
a subset of which are the first class constraints <f$ . These first class constraints 
are used in the formation of the quotient space 

M i 

p i = 77T • 

The projectability condition for Hi (the pullback of H c to Mi) to Vi may be expressed as 
the condition {H c , (f)^} |mi = 0. If this condition holds, we have found the final physical 

(2) 

space. If it doesn't, there are new, secondary constraints 0^ , and these constraints 
along with the initial set of primary constraints 0jL are used to define a constraint 
surface M 2 . Among the set of constraints used to define M 2 are first class constraints, 
including a subset of the first class constraints associated with Mi, which we denote 
0^2, and a subset of the secondary constraints, which we denote 4$. These first class 
constraints are used in the formulation of the quotient space 

2 '~ ({-, </>£},{-, 02}) ' 

Again we must require projectability of the Hamiltonian; eventually, the final phase 
space is 

M f 

,(!)-. r "7(27-, , a-., 



({- 4>V f }> {-, fa,}, ■■■,{-, Mi }) 

where (f>^) are the final first class n-ary constraints, all of which are taken in effective 
form. Vf is endowed with a symplectic form which is the projection of the form ujf 
in Mf, which is the final constraint surface. The dimension of Vf is 2N — M — Pf, 
where iV is the dimension of the initial configuration space, M is the total number 
of constraints, and Pf is the number of final first class constraints. Observe that we 
end up with the standard counting of degrees of freedom for constrained dynamical 
systems: First class constraints eliminate two degrees of freedom each, whereas second 
class constraints eliminate only one each. The final result is an even number because 
the number of second class constraints is even. 

In order to use the technique we've presented, the constraints are presumed to 
be effective (for example, see equation (|2.4j )) — if ineffective constraints occur, they 
can always be made effective for use with this technique; in that sense, the technique is 



Reduced phase space 



11 



actually geometrical. One might ask whether such modification of ineffective constraints 
can cause problems. It turns out that if ineffective constraints occur, then their presence 
may modify the gauge fixing procedure used in conjunction with the original Dirac 
method in such a way that the counting of degrees of freedom differs from the one 
presented above. In the next section we discuss a simple example that shows the 
difference between Dirac's original treatment, supplemented by gauge fixing, and the 
quotienting method we've outlined here, which corresponds to Dirac's extended method. 

Dirac's extended method, which is equivalent to the one we've presented here, 
produces a final phase space which is always even dimensional. Dirac's original 
procedure, supplemented by gauge fixing, has the superiority of being equivalent to 
the Lagrangian variational principle. Therefore, in spite of the fact that this latter 
method may result in a system with an odd number of degrees of freedom (as in the 
example in the following section), it is to be preferred for classical models. 

5. Example 

Consider the Lagrangian 

£ = ii 2 + ^> 2 , (5.D 

where 2^0. The Noether gauge transformations are 

5x — , 5y = — , Sz = e , 

z 

where e is an arbitrary function. 

First, we analyze this system from a Lagrangian point of view. The equations of 
motion are 

x = , y = . (5.2) 

The z variable is completely arbitrary and is pure gauge. These equations define a system 
with three degrees of freedom in tangent space, parameterized by x(0),x(0),y(0). Notice 
that the gauge transformation 5y vanishes on shell, so y is a weakly gauge invariant 
quantity. 

Let us now analyze this system using Dirac's method. The Dirac Hamiltonian is 

Hd = -pl + -zp 2 y - Xp z , (5.3) 

where A is the Lagrange multiplier for the primary constraint p z = 0. The stabilization 
of p z = gives the ineffective constraint p 2 y = 0, and the algorithm stops here. The 
gauge generator (Batlle et al 1989, Pons et al 1997) is 

G = ep z + e -p\ , (5.4) 

with e an arbitrary function of time. 

The gauge fixing procedure, Pons and Shepley (1995), has in general two steps. The 
first is to fix the dynamics, and the second is to eliminate redundant initial conditions. 
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Here, to fix the dynamics we can introduce the general gauge- fixing z — f(t) = for / 
arbitrary. Stability of this condition under the gauge transformations sets e = 0. Since 
the coefficient of e in G is ineffective, it does not change the dynamical trajectories, and 
so the gauge fixing is complete. Notice that this violates the standard lore, for we have 
two first class constraints, p z = and p y = but only one gauge fixing. This totals three 
constraints that reduce the original six degrees of freedom to three: x(0),p x (0), y(0), the 
same as in the Lagrangian picture. 

Instead, if we apply the method of quotienting out the kernel of the presymplectic 
form, we get final reduced phase space 
_ M 2 

f ({- P*M- Py}) ' 
where M 2 is the surface in phase space defined by p z = 0,p y = 0. We have uj 2 = dxAdp x 
and H 2 = \p x . The dimension of Vj is 2. 

This result, which is different from that using Dirac's method, matches the one 
obtained with the extended Dirac's Hamiltonian, where all final first class constraints 
(in effective form) are added with Lagrange multipliers to the canonical Hamiltonian. 
Dirac's conjecture was that the original Dirac theory and the extended one were 
equivalent. We conclude that when Dirac's conjecture holds, the method of quotienting 
out the kernel is equivalent to Dirac's, whereas if Dirac's conjecture fails, it is equivalent 
to the extended Dirac's formalism. 

6. Conclusions 

In summary, we have the following. 

1) We have obtained a local basis for /C = Ker(cjL) in configuration-velocity space 
for any gauge theory. This is new and allows for trivial verifications of the properties of 
K given in the literature. To get these results it has been particularly useful to rely on 
Hamiltonian methods. 

2) We have obtained as the final reduced phase space the quotient of the final 
Dirac's constraint surface in canonical formalism by the integral surface generated by 
the final first class constraints in effective form. We find the constraint surface (Mf 
in equation ( f4.4| )) through a projectability requirement on the Lagrangian energy (or 
equivalently, on the Hamiltonian) rather than through imposing tangency conditions on 
the Hamiltonian flows. Let us emphasize this point: We do not talk of stabilization of 
constraints but rather projectability of structures which are required to formulate the 
dynamics in a reduced physical phase space. 

3) We have compared our results with Dirac's procedure. An agreement exists in all 
the cases when no ineffective Hamiltonian constraints appear in the formalism. If there 
are ineffective constraints whose effect iviz at ion is first class, then our results disagree 
with Dirac's, and it turns out that the quotienting algorithm agrees with the extended 
Dirac formalism. When there are disagreements, the origin is in the structure of the 
gauge generators. Sometimes the gauge generators contain pieces that are ineffective 
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constraints, and they, contrary to the usual case, do not call for any gauge fixing. 
Essentially, the variables that are canonically conjugate to these first class ineffective 
constraints are weakly (on shell) gauge invariant. The quotienting reduction method, 
as well as Dirac's extended formulation, eliminates these variables and yields a phase 
space whose variables are strictly (on and off shell) gauge invariant. This is the difference 
with Dirac's original method, supplemented with gauge fixing, which is able to retain 
the weakly gauge invariant quantities. For this reason we feel that this latter technique 
is superior to the quotienting algorithm in these circumstances — at least for classical 
models. 

4) We have produced a simple example that illustrates the failure of Dirac's 
conjecture in the presence of ineffective constraints. This example also shows that, 
in Dirac's analysis, it is possible to have Hamiltonian formulations with an odd number 
of physical degrees of freedom. We must remark that in Dirac's approach (supplemented 
with gauge fixing) it is not always true that every first class constraint eliminates two 
degrees of freedom: This does not happen if there are first class constraints that appear 
in the stabilization algorithm in ineffective form. 

5) It is worth mentioning that other reduction techniques, specifically the Faddeev 
and Jackiw method, may also fail to reproduce Dirac's theory (Garcia and Pons 1998) 
if the formalism contains ineffective constraints. 

6) Of course, one should not forget quantum mechanics. The canonical approach 
to quantum mechanics involves a (nonsingular) symplectic form, Isham (1984). In this 
method, it is therefore required that phase space be even-dimensional. This argument 
would tend to favor the quotienting algorithm. However, it may be that other approaches 
to quantum mechanics, possibly the path integration approach, do not need such a 
requirement. And in any case, it is not strictly necessary that a model which is 
acceptable as a classical model be quantizable. It is for these reasons that we say that 
an approach to Hamiltonian dynamics which results in a phase-space picture equivalent 
to the tangent space picture — the original Dirac method supplemented with gauge 
fixing — is favored for classical models. 
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